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Abstract 

We consider the Born-Oppenheimer problem near conical intersection in two 
dimensions. For energies close to the crossing energy we describe the wave 
function near an isotropic crossing and show that it is related to generalized 
hypergeometric functions o-^^s- This function is to a conical intersection what 
the Airy function is to a classical turning point. As an application we calculate 
the anomalous Zeeman shift of vibrational levels near a crossing. 

PACS numbers: 31.15.Gy, 33.55.Be, 33.20.-t 



1 



Typeset using REVTgX 



Introduction. The Born-Oppenheimer (BO) problem is concerned witli tlie analysis 
of Schrodinger type operators where the small electron to nucleon mass ratio, plays the 
role of the semiclassical parameter. The theory identifies distinct energy scales: The 

electronic scale which, in atomic units, is of order one and the scale of nuclear vibrations 
which is of order (l/M)^/^ in these units. M is the nucleus to electron mass ratio. The 
identification of the electrons as the fast degrees of freedom is central to the theory. 

The clean splitting between fast and slow degrees of freedom fails near eigenvalue crossing 
of the electronic Hamiltonian where there is strong mixing between electronic and vibrational 
modes. This lies at the boundary of the conventional BO theory. Since the coupling between 
different electronic energy surfaces becomes infinite near a crossing, the nuclear wave function 
does not reduce to a solution of a scalar (second order) Schrodinger equation. We describe 
the (double surface) nuclear wave function near an isotropic conical crossing, for energies 
close to the energy of the crossing. 

The strong mixing of the electronic and nuclear degrees of freedom near crossing leads 
an anomalous Zeeman effect. To describe the anomaly recall that the Zeeman splitting 
in molecules is reduced compared to the Zeeman splitting of atoms. It is convenient to 
parameterize the reduction by a parameter 7 so that the Zeeman splitting is of the form 
(and order) M~^B with B the external magnetic field. The low lying vibrational levels have 
large reduction, 7 = 1. This is what one expects from nuclei whose magnetic moments are 
by factor M smaller than the Bohr magneton. Levels near the crossing energy can have a 
small reduction expressed by the fact that 7 < 1. For the isotropic situation we calculate 
7 = 1/6, so that the Zeeman shift is anomalously large, by a factor of about 2000, than 
the normal Zeeman splitting of molecular levels. More precisely, we find that the Zeeman 
splitting near isotropic crossing is 

AE{m)^-^^-^B. (1) 

The sign means equality in the limit M — > 00. B is proportional to the magnetic field, with 
a coefficient dependent only on the electronic wave functions at the crossing, m, a half odd 
integer, is the azimuthal quantum number, is an electronic time scale, see Eq. (plQ]) below. 
g{m) is a universal dimensionless factor which is determined by the nuclear wave function 
near the crossing, see Eq. ([T8|) below. As we shall see g{m) = —g{—m) and numerical 
estimates of Eq. (|l^) give 

g = 0.961, g (^^) = 0.543, g (^^) = 0.396. (2) 

g{m) is a molecular analog of the Lande g factor in atoms: So, while Lange g factor describes 
the Zeeman shift due to the mixing of spin and orbital degrees of freedom, g{m) does it for 
the nuclear and electronic ones. 

One can formulate the BO problem in the following way [§]: 

Hbo = -j^A, + H,{x), (3) 
where Hg{x), the electronic Hamiltonian, depends parameterically on the nuclear coordinates 



X. When time reversal is not broken, Hf,[x) is a real symmetric matrix ||T0[- The Wigner von 
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Neumann crossing rule says that He{x) has generically a crossing point for two modes 
of vibrations, x G R^. 

Here we shall consider the simple scenario where He(x) is a 2 x 2 matrix of and x G R^. 
This means that we shall treat only the restriction of the electronic Hamiltonian to the two 
dimensional subspace spanned by the two degenerate eigenstates at the crossing point. We 
shall assume that He{x) has a single crossing point at a; = 0, and set the crossing energy 
at 0. We further assume that the crossing is conic, that He{x) is isotropic about the origin, 
and that the x dependence of Hf.{x) is smooth near the origin. This is a common model [0. 

We first recall why the standard BO theory fails near a crossing. When H(,{x) is symmet- 
ric it can be diagonalized by an ort/iog'ona/ transformation R{x). In the 2x2 case, and when 
H(.{x) is non-degenerate, R{x) is uniquely determined, up to an overall sign, by requiring 
det R{x) = 1. Hence, away from crossing, Hbo is unitarily equivalent to 

where Ei2{x) are the two eigenvalues of Hf.{x). The vector potential A{x) = 
—iR^ {x)'V xR{x) is purely off-diagonal. For linear crossing, R{x) —R{x) as x surrounds 
the origin This forces a l/\x\ singularity of the vector potential for small x. Far from 
the origin, to leading order in 1/M, the two components of the wave function decouple and 
are given by 

^boA^)^iJci,j{^)R{x)(^^';'), J = 1,2, (5) 



where ipcij is a semiclassical solution of the Schrodinger equation with potential Ej{x) in 
the cut plane with antiperiodic boundary conditions @J^. This decoupling holds provided 
>> M~^/^ [^. We shall henceforth denote e = M^^^^, and refer to the region |x| >> e 



\x 



as "far from the crossing". In contrast, the divergence of the off diagonal part of A near 
the crossing prevents such a decoupling near the crossing. Our aim is to analyze the BO 
theory, to leading order in e, near the crossing of He- The reason one can still hope to say 
something useful near the crossing is that the asymptotic form of He{x) near the crossing, 
i.e. for << 1, is universal 

He{x) = xiai + X2a'i + 0{x'^) (6) 

where a are the Pauli matrices. We shall refer to |a;| << 1 as "close to the crossing". Notice 
that our notions of close and far from the crossing have a nonempty intersection. This 
enables us to match the solution close to the crossing with the standard, decoupled, BO 
solution, Eq. (|^). 

The zero energy wave function close to the crossing. We study the wave functions 
near the crossing for energies close to the crossing energy. Everything to be said from now 
on is true in the limit of M — >■ cx), to leading order in negative powers of M. 

We shall assume that zero is an eigenvalue of (|^), since there is always an eigenvalue 
that is close enough to zero [Q. It turns out to be convenient first to unitary-transform 
with e~*4'^^ This will replace as in (D) by cr2- In this representation, close to the crossing 
the zero energy wave function satisfies approximately the differential equation 
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{-V^ + eiCTi +60^21^(6,6) = (7) 

\1/ stands for a two component column matrix and 6 = M ' scaled variable. 

The operator J3 = L3 + |cr3 = —i^id^^ + i6f^5i + |'^3 commutes with the operator on 
the l.h.s. of (^. It does not have the meaning of total angular momentum, since the Pauli 
matrices do not represent spin. We thus consider solutions of the differential equation (^ 
which are eigenf unctions of J3 with an eigenvalue m, namely: 



Amd 



ipi{p;m)e 



(8) 



p, 6 are the polar coordinates related to 6 ^ must be half odd integer, for the wave function 
to be single valued. Separating variables, the radial equation obtained from (|^), in the m-th 
sector, takes the form: 



dfP' p dp 



Id m? 
+ 



P 




(9) 



Of the four linearly independent solutions of Eq. (^, basically due to boundary conditions, 
only one linear combination is fit to represent a wave function. We denote it by J-'c{m, p) 
and its components by ipidm; p) and ip2c{m;p). J-'c{m, p) is to a crossing point what the 
Airy function is to a classical turning point |]I2|,|13|: It interpolates between the near region 



where the wave function is intrinsically a two component spinor, and the far region where 
the wave function is highly oscillatory and given by Eq. (^. J-'djn, p) has a closed expression 
in terms of the generalized hypergeometric functions of the kind 0-^3- It has the asymptotic 
(for p >> 1) form 



^c{m, p) 



P 



1 /2 fm 1 



(10) 



Solving (pp asymptotically near the origin gives 15 



V9i(p; m) 
V'2(p; m) 



,m-l/2 

6+p'"+i/2 



a^p"" + a_p 



-m+1/2 



(11) 



Solving (P) asymptotically at infinity we obtain 



y?i(p;m) 
V52(p; m) 



1 



P 



,3/4 



C cos [z + 



-1 



;i2) 



The four dimensional family of solutions can be parameterized by either a+,a_,6+,6_ or 
Aj^, A_, C, (f). Requiring the solution to be bounded at the origin and at infinity means (for 
positive m) that A^ = 0, a_ = 0, h_ = 0. Imposing three homogeneous conditions on a four 
dimensional linear space leaves us with one dimensional subspace, i. e. a certain function 
times an arbitrary constant. This is the celebrated J-'c{p,m). 

While the reason we require J-'c{p,m) to be bounded at the origin is obvious, the reason 
we require it to be bounded at infinity is a bit subtle, since (^ is meaningful only close to 
the crossing. However, "close to the crossing", p « M^^^, extends farther and farther in 
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terms of p as M gets large. at ([T2|) should be exponentially small, and can set to zero 
to leading order. 

The solutions which are regular at the origin can be obtained from the fourth order 
differential equations, one for each component (pj, obtained from (|^). These are related to 
the differential equation that defines the generalized hypergeometric functions o-^s- The two 



linearly independent solutions that are regular at the origin are |20 

J^Jp- m) = I m^o^^l, 3' 2 ^ T' 6 ^ I (13) 



where o-^sl; 6, c; p) are generalized hypergeometric functions [p!6| , p!7| . The linear combina- 
tion 

J^c{p; m) = Ai{m)J^i{p; m) + A2{m)J^2{,p; m) (14) 
is bounded at infinity provided [Q: 

3r(i + f )r(i - ti)i)r(i + ^ ' 



We have set the global coefficient in front of (14) such that ([l^) will be true. Formulae 



( p!3| - p!5| ) are correct only for positive m-s. Their negative counterparts can be obtained by 



interchanging the lower and upper components, as can be seen from (^. 

The anomalous Zeeman shift. Let us now turn to the Zeeman shift. The magnetic 
field breaks time reversal symmetry, which means that it adds an the imaginary term to 
He{x), in the representation where He{x) is real. Generically the magnetic field will remove 
the conical intersection of He{x) at a; = and create a gap between the two energy sheets. 
The gap will be proportional to the magnetic field in atomic units, and the coefficient will 
be in general of order one. We can therefore introduce the magnetic field into our model 
by adding the term Ba2 to He{x) at @. There is no harm in taking B to be independent 
of X, since only the value of B in the origin will be of importance to us. B is therefore a 
constant, proportional to the magnetic field in atomic units. We shall not minimally couple 
the magnetic field directly to the vibrations for this turns out to be a weaker effect, of order 
M~^, while the shift mediated by the electrons, in the rotationally invariant case, is M~^/^, 
as we shall see. 

Let us consider the case where the rest of He{x), namely the O(x^) part in (^, does 
not break the full rotational symmetry of its linear part, so that m is still a good quantum 
number. 

In the rotational invariant case, the model has, for B = 0, a. two-fold degeneracy: The 
states m and —m are degenerate. The magnetic field B breaks this degeneracy. The splitting 
is twice the Zeeman shift in the energy for the two states ±m move in opposite directions. 

Equipped with approximants to the wave function near and far from the crossing, with 
a degenerate perturbation theory one can calculate, to leading order in 1/M, the Zeeman 
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splitting and obtain (|l|). We describe this calculation in details elsewhere |2^. Here we 
would only like to sketch the derivation of the power in (|l]). 

Far from the crossing, one can neglect the vector potential in (^. The WKB approximant 
to the radial part of ^cz,2 is 

N 2 

where we have employed the linearity of the energy surfaces near the crossing. It is a general 
property of the WKB approximation, that the normalization coefficient N is independent 
of M, to leading order in 1/M. From ([T0|) and (|16D one sees that the interpolation of the 



BO radial wave function towards the crossing is 

^(r; m) ^ NM^^^J^dm; M^''^r). (17) 

B removes the degeneracy of the electronic levels at x = 0. The gap created there due to B 
is equal to 2B. Intuitively, the Zeeman shift of a vibrational level will be proportional to the 
amount of probability density in the vicinity of the crossing times B. By "vicinity" we mean 
a neighbourhood of order M~^/^ of the origin, the area of which is of order M~^/^. The 
density associated with the wave function is large there and by (|1^) is proportional to M^/^. 
Hence the weight is proportional to (M)~^/^ which gives the power of the Zeeman splitting 
in Eq. (Q). The coefficient of proportion will include an integral over the components of JF, 
which gives g{m) 



g{m) = pdp{^^^{p;m) -ip2^{p;m)). (18) 

JO 

N gives the factor l/Tg, where Tg is 



T, = / dr (19) 



The integration is carried out between the two turning points of Ei, the negative energy 
sheet (§). Eq. (|^) is proportional to the time in takes a particle with a unit (electronic) 
mass to travel classically across the potential. It is independent of the nuclear mass, and 
T^^ has the order of magnitude of electronic energies. From the invariance under T of 
Hbo{B = 0) it follows that g{m) = —g{—m). 

One motivation for this work was an attempt to gain some understanding of the differ- 
ent status of crossing in theory and experiment. Theory puts crossing and avoided cross- 
ing in distinct baskets: conic crossings come with fractional azimuthal quantum numbers 
while avoided crossings come with integral quantum numbers. In contrast, measurements 
of molecular spectra normally can not tell a crossing from near avoided crossing. Only with 



precision measurements and precise quantum mechanical calculations ||19[ can one tell 
when molecular spectra favor an interpretation in terms of crossing and half integral quan- 
tum numbers or avoided crossing with integral quantum numbers. Zeeman splitting appears 
to be a useful tool to study crossing. The anomaly of crossing is characterized by a fractional 
power 7 of the molecular reduction of the Zeeman splitting, M"'''. A system of choice for 
studying crossing is molecular trimers Since trimers are not rotationally invariant 
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Eq. (HD does not apply and we can not conclude that near crossing 7 = 1/6 for trimers. 
However, it is natural to expect that the qualitative features of our results carry over also to 
the non-isotropic case where crossing will manifest itself by an anomalously large Zeeman 
splitting and a fractional 7. It is an interesting challenge to calculate or measure the value 
of 7 for (other) molecular crossings, and trimers in particular. 
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FIGURES 




FIG. 1. The components of !Fci'm'', p) for m = ^ and their asymptotic form given by (|lO|), which 
is also the BO wave function approximated by WKB. 
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